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Abstract 

This  paper  develops  a  macroscopic  polarization  switching  model  which  characterizes  the  ferroelastic 
switching  mechanisms  inherent  to  lead  zirconate-titanate  (PZT)  in  a  manner  suitable  for  subsequent 
transducer  and  control  design.  We  construct  Helmholtz  and  Gibbs  energy  relations  at  the  lattice  level 
which  quantify  the  internal  and  electrostatic  energy  associated  with  90°  and  180°  dipole  orientations. 
Equilibrium  relations  appropriate  for  homogeneous  materials  in  the  absence  or  presence  of  thermal 
relaxation  are  respectively  determined  by  minimizing  the  Gibbs  energy  or  balancing  the  Gibbs  and  rel¬ 
ative  thermal  energies  using  Boltzmann  principles.  Macroscopic  models  suitable  for  nonhomogeneous, 
polycrystalline  compounds  are  constructed  through  stochastic  homogenization  techniques.  Attributes 
and  limitations  of  the  model  are  illustrated  through  comparison  with  experimental  PLZT  data. 

1  Introduction 

This  model,  extending  the  work  done  in  [1],  combines  energy  principles  at  the  lattice  level,  theory 
of  thermally  activated  processes,  and  stochastic  homogenization  techniques  to  characterize  hysteresis 
due  to  ferroelectric  and  ferroelastic  switching  in  a  manner  which  facilitates  material  characterization, 
transducer  design,  and  model-based  control  design.  Helmholtz  and  Gibbs  energy  relations  are  con¬ 
structed  at  the  lattice  level  to  quantify  the  internal  and  electrostatic  energy  associated  with  90°  and 
180°  dipole  orientations.  For  regimes  in  which  thermal  activation  is  significant,  the  Gibbs  and  relative 
thermal  energies  are  balanced  through  Boltzmann  theory  to  provide  equilibrium  relations  quantifying 
local  strains  and  polarizations  as  a  function  of  input  stresses  and  fields.  In  the  limit  of  negligible 
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thermal  activation,  the  local  relations  reduce  to  minima  of  the  Gibbs  energy  and  enforcement  of  these 
minimization  criteria  significantly  improves  the  efficiency  of  the  model  and  inversion  process  employed 
for  linear  control  design.  Parameters  such  as  the  local  coercive  and  interaction  fields  are  considered 
to  be  manifestations  of  underlying  distributions  to  incorporate  the  effects  of  polycrystallinity  and  ma¬ 
terial  nonhomogeneities.  Homogenization  in  this  manner  yields  low-order  models  posed  in  terms  of 
effective  parameters  that  are  efficient  to  implement  in  optimization  or  control  algorithms. 

Physical  mechanisms  associated  with  ferroelectric  and  ferroelastic  switching  are  summarized  in 
Section  2  to  motivate  90°  and  180°  switching  mechanisms  that  need  to  be  incorporated  in  energy 
relations.  In  Section  3,  we  summarize  the  construction  of  a  piecewise  defined  Gibbs  energy  functional. 
Minimization  of  this  functional  quantifies  homogeneous,  single  crystal  behavior  in  the  absence  of 
relaxation  processes.  In  Section  4,  stochastic  homogenization  techniques  are  employed  to  construct  a 
macroscopic  model  for  polycrystalline,  nonhonrogeneous  materials  and  in  Section  5,  properties  of  the 
model  are  illustrated  through  comparison  and  prediction  of  experimental  PLZT  data. 

2  Ferroelectric  and  Ferroelastic  Switching  Mechanisms 

Lead  zirconate-titanate  is  comprised  of  PbZr^-iOa  (lead  zirconate)  and  PbTi^Os  (lead  titanate)  where 
x  is  chosen  to  optimize  electromechanical  coupling.  For  temperatures  above  the  Curie  point  Tc,  the 
structures  of  PbTiOa  and  PbZr03  are  cubic  whereas  for  T  <  Tc,  the  structure  of  PbTiOs  is  tetragonal 
and  PbZr03  is  orthorhombic  [6].  The  switching  mechanisms  are  illustrated  in  the  context  of  the 
paraelectric  cubic  and  the  ferroelectric  tetragonal  structure  of  lead  titanate  and  note  the  analogous 
behavior  is  observed  for  the  orthorhombic  structure  of  lead  zirconate. 

The  application  of  an  electric  field  E  that  is  larger  in  magnitude  than  the  coercive  field  Ec  induces 
ferroelectric  switching.  For  PbTi03,  this  causes  the  central  Ti+4  ion  to  relocate  to  a  new  equilibrium 
position,  resulting  in  a  180°  change  in  polarization  that  is  parallel  to  the  applied  field  as  depicted  in 
Figure  1(a).  Ferroelastic  switching  is  caused  by  the  application  of  a  stress  a  that  is  larger  in  magnitude 
than  the  coercive  stress  ac  producing  a  90°  change  in  polarization  that  is  perpendicular  to  the  applied 
stress  as  illustrated  in  Figure  1(b).  The  ferroelectric  and  ferroelastic  switching  mechanisms  cause 
a  hysteretic  relationship  between  input  fields  E  and  a  and  output  polarization  P  and  strains  e  as 
illustrated  in  Figure  2. 

At  point  A  in  Figures  2(a)  and  2(b),  the  electric  field  is  sufficiently  strong  so  that  all  the  dipoles  form 
one  domain  that  is  aligned  in  the  direction  of  the  applied  field.  As  the  field  is  decreased,  it  approaches 
the  negative  coercive  field  in  the  region  around  point  B  where  180°  switching  commences.  Additionally, 
this  often  includes  90°  switching  as  indicated  by  the  presence  of  a  negative  strain  at  point  B  as  depicted 
in  Figure  2(b).  At  point  C,  all  the  dipoles  have  switched  and  again  form  one  domain  that  is  aligned  in 
the  direction  of  the  electric  field.  At  point  C,  the  polarization  is  opposite  to  that  at  point  A  whereas 
the  strains  have  the  same  value.  As  the  field  is  increased,  90°  switching  occurs  at  point  D  and  rapidly 


(a) 


(b) 


Figure  1:  (a)  Ferroelectric  180°  switch  in  spontaneous  polarization  Pq  induced  by  an  applied  electric 
field,  and  (b)  ferroelastic  90°  switch  induced  by  an  applied  stress. 
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Figure  2:  (a)  Hysteretic  field-polarization  relation  for  bulk  PLZT,  (b)  field-strain  behavior  of  PLZT, 
(c)  stress-polarization  relation  for  PLZT,  and  (d)  stress-strain  behavior  of  PLZT. 


continues  back  to  point  A  where  the  full  180°  switch  has  occurred  and  the  dipoles  are  again  aligned 
with  the  applied  electric  field.  Figures  2(c)  and  2(d)  illustrate  the  relationship  between  an  externally 
applied  stress  and  the  polarization  and  strains.  At  point  F,  the  dipoles  are  aligned  in  the  positive 
3-direction  and  the  material  acts  as  one  domain.  As  the  compressive  stress  is  increased  in  magnitude, 
it  approaches  the  coercive  stress.  In  the  region  around  point  G,  90°  switching  occurs  and  the  dipoles 
begin  to  align  perpendicular  to  the  direction  of  the  applied  stress.  This  is  indicated  by  the  presence 
of  a  negative  strain  at  point  G  in  Figure  2(d).  As  the  stress  is  reduced  in  magnitude,  the  material 
remains  poled  perpendicular  to  the  applied  stress  resulting  in  a  decrease  in  the  polarization  in  the 
3-direction. 


3  1-D  Ferroelastic  Switching  Model 

To  facilitate  model  implementation,  improve  computational  efficiency  and  increase  the  model  flexibil¬ 
ity,  we  employ  piecewise  quadratic  energy  functionals  which  are  approximations  of  energy  functionals 
based  on  Landau-  Devonshire  principles.  The  model  coefficients  directly  relate  to  measured  material 
properties  allowing  a  means  for  parameter  identification  and  estimation.  The  low-order  piecewise 
polynomial  models  also  facilitate  transducer  and  control  design  and  incorporates  both  90°  ferroelec¬ 
tric  and  ferroelastic  polarization  switching  mechanisms.  We  note  that  aspects  of  the  functionals  are 
similar  to  those  employed  for  SMA  undergoing  austinite  -  martinsite  phase  transformations  and  the 
reader  is  referred  to  [3,  6,  7]  for  details  illustrating  properties  of  the  SMA  relations. 


3.1  Helmholtz  and  Gibbs  Energy  Relations 


We  consider  electric  field  and  stress  inputs  (E,a)  and  polarization  and  strain  outputs  (P,  e).  To 
account  for  the  180°  and  90°  polarization  switching,  the  polarization  has  three  allowed  dipole  states 
P_,  P+  and  -Pgo-  We  define  the  polarization  component  of  the  Helmholtz  energy  to  be 
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i 

VI 
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where 


p  _  V  {Pi  -  Pr)  P901  -  V2P901P1  _  Pi  ~  Pr 

m~  r,  (Pj  -  PR)  -  rnPm  ’  Vl  ~  V  Pi  -  Pj  (1) 

P=\{Pi~  Pr)2  ~  |  (Pi  -  Pm)2  ,  A=|  (P901  ~  Pm)2  +  P  ~  f  (^90/)2  • 

The  electromechanical  coupling  component  is  given  by 

'ipesiP,  e)  =  -aeP  -  (2) 

where  a  is  the  piezoelectric  coupling  coefficient  and  q  is  the  electrostrictive  coupling  coefficient.  The 
elastic  energy  is 

Al(s)  = 

where  Y  is  the  Young’s  modulus.  The  total  Helmholtz  energy  is  then  given  by 

ip{P,  e)  =  'ipp(P)  +  ipes(P,  e)  +  V’ez(e)- 

Balancing  the  internal  energy  ip  and  the  externally  applied  energy  yields  the  Gibbs  energy 

G(E,  P,  a,  s)  =  ip(P,  e)  —  EP  —  as.  (3) 

The  necessary  condition  ^  =  0  yields 

s=Y~1  ( a  +  aP  +  qP 2)  . 

In  a  manner  analogous  to  that  employed  in  [4]  for  a  stress-free  state,  the  strain  is  substituted  into  the 
Gibbs  energy  (3)  resulting  in  a  Gibbs  energy  of  the  form 

G(E,  P,  a)  =  $(P,  a)  —  EP  —  ae.  (4) 


3.2  Switching  in  the  Absence  of  Thermal  Activation 


Neglecting  thermal  effects,  the  local  polarization  ( P )  for  a  single  crystal  with  an  applied  stress  can  be 
determined  by  solving  the  necessary  condition  =  0,  or  by  minimizing  the  Gibbs  energy 

(P)  =  arg  nrin  G(E,  P,  a)  (5) 


where  the  Gibbs  energy  is  defined  by  equation  (4). 

To  calculate  the  value  of  the  coercive  field,  we  apply  the  necessary  condition  §p  =  0  and  solve  for 
the  value  of  the  electric  field  that  results  in  P  =  Pj  to  obtain 


Ec(a) 


‘IqPi-a  2  q2Pf 


3  aqPj 

Y  + 


a2  Pi 
Y 


+  V  (Pr  -  Pi) . 


(6) 


It  is  clear  from  (6)  that  the  coercive  field  is  stress  dependent.  The  coercive  stress  ac  can  be  computed 
by  setting  Ec  =  0  in  (6).  This  results  in 


<Tc  =  -  1  \jfY (Pr  ~  Pi)  +  2q2(Pjf  -  ZaqiPj)2  +  a2(P/)]  . 


(7) 


The  electric  field  required  to  eliminate  the  90°  minima  in  the  Gibbs  energy  is  denoted  by  E®° .  An 
explicit  expression  for  the  value  of  E®°  is  found  to  be 


El\a)  =  -(a) 


a  +  2 qP90  2 q2Pg0  +  3 aqPg0  +  (a2  -  rj2Y)P9 0 


Y 


Y 


(8) 


A  direct  180°  ferroelectric  polarization  switch  (i.e.,  P±  — >  P^)  can  occur  if  E^°  <  Ec.  A  90°  ferro¬ 
electric  polarization  switch  (i.e.,  P±  — >  P99)  can  occur  if  E^°  >  Ec  and  a  90°  ferroelastic  polarization 
switch  will  occur  when  a  >  ac. 
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3.3  Switching  in  the  Presence  of  Thermal  Activation 

For  the  switching  criteria  summarized  in  Section  3.2,  dipole  switching  occurred  only  when  the  magni¬ 
tude  of  the  applied  electric  field  or  stress  exceeded  that  of  the  coercive  field  or  stress  value  of  the  single 
crystal.  However,  mechanisms  such  as  excitation  from  thermal  effects  can  induce  switching  before 
local  minima  are  eliminated.  To  model  this  phenomena,  we  incorporate  the  effects  of  thermal  activa¬ 
tion.  We  continue  to  assume  that  the  material  is  homogeneous  and  that  each  dipole  in  this  collection 
has  the  same  energy  landscape  and  are  subject  to  the  same  switching  mechanisms.  This  allows  us  to 
model  the  switching  that  occurs  by  evolving  the  fraction  of  dipoles  in  each  allowed  orientation  by  a 
population  model  that  is  similar  to  what  is  derived  in  [8]. 

In  the  manner  detailed  in  [6] ,  thermal  activation  processes  are  incorporated  through  the  use  of  the 
Boltzmann  relation 

fi(G(P))  =  Ce~G{p)v/kT.  (9) 

Here  V  denotes  a  representative  control  volume,  k  is  Boltzmann’s  constant  and  T  is  the  temperature. 
The  fraction  of  dipoles  in  each  allowed  orientation  is  given  by  x-,xgo  and  x+.  Conservation  of  the 
number  of  dipoles  yields 

X-+Xgo  +  X+  =  l.  (10) 

The  local  polarization  (P)  is  given  by 


(P)  =  x-  (P-)  +  Xgg  ( Pg0 )  +  x+  (P+) 


(11) 


where  (P-) ,  (P90)  and  (P+)  are  the  expected  polarization  values  associated  with  each  allowed  dipole 
orientation.  As  illustrated  in  [6,  8] ,  these  values  are  found  by  integrating  the  product  of  the  polarization 
P  and  the  Boltzmann  probability  density  p(G(P))  over  the  allowed  polarization  states.  This  simplifies 
to  the  relations 


(P-) 


f~Pl  pe~7G(P)dp 

J  —  OO 

r-Pl  e-^G(P)dP  ’ 

J  —  OO 


(P9o) 


I- 


p90!  pe-lG{P)dp 

*901 


rp9oi  e-lG(p)dP 
J  —  Tqn  T 


(P+) 


joo  pe-^G{P)dp 

e~^G(P)dP 


where  7  =  .  For  computational  efficiency,  we  can  approximate  these  values  by  use  of  the  necessary 

condition  pp  =  0  which  yields 


(P-)  =  --Pr,  (Pg  0)  =  -,  (P+)  =  -+PR. 
n  rj2  T] 


In  this  case,  (P_) ,  (Pgg)  and  (P+)  are  the  location  of  the  minima  of  the  Gibbs  energy. 

The  time  evolution  of  the  dipole  phase  fractions  are  governed  by  the  first  order  ODE  system 


(  X-  \  (  -P-  P90—  0  \  /  x-  \ 

j  Xgg  J  =  p-  -(P90-+P90+)  P+  Xgg 

\  X+  )  V  0  P90+  -P+  J  \  X+  J 


The  system  results  from  the  assumption  that  transitions  between  the  three  allowed  states  occur  only 
to  the  nearest  neighbor  and  can  be  reduced  in  dimension  by  applying  the  conservation  relation  (10). 
The  likelihood  to  switch  out  of  the  P-  orientation  into  the  P90  orientation  is  denoted  by  p-  and  the 
notation  for  the  remaining  likelihoods  is  summarized  in  Table  1. 

The  likelihood  p-  is  calculated  by 

1  e-7  G(-Pi) 

P ~  ~  r  [~pi  e-'rG(P)dP 

J —00 
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State  Transition 

notation 

P_  - 

4  -P90 

P- 

P9  0  - 

->  R_ 

P90— 

Pm  ~ 

-  I\ 

P90+ 

P+  - 

4  -P90 

P+ 

Table  1:  1-D  transition  likelihoods. 


where  the  relaxation  time  r  is  the  reciprocal  of  the  frequency  at  which  dipoles  attempt  a  switch.  The 
likelihood  of  switching  out  of  the  P90  orientation  into  the  P  orientation  is  specified  by  pm-  and  is 
calculated  by 

1  g— 7G(— -P90) 

P9°~  ~~  r  jP^oe-iG(P)dp- 

The  likelihoods  P90+  and  are  obtained  in  a  similar  manner.  The  likelihoods  can  also  be  evaluated 
in  terms  of  error  functions  resulting  in, 

1  g-7  G(P/)  !  e-a(Pjr+6/: 2)2 

P+  ~  t  e~^pUP  ~  T\^/l.exic(^(PI  +  b/2)) 


where  a 


rn 

2 


and  b  =  —2 


(Pr  ~  E ) 


V 


The  likelihoods  P90+ , P90-  and  p-  follow  similarly. 


4  Macroscopic  Polarization  Model 

Nonuniformities  in  the  lattice  structure  due  to  polycrystallinity,  material  nonhomogeneities  and  vari¬ 
ations  across  grain  boundaries  produce  a  distribution  of  Helmholtz  and  Gibbs  energy  profiles  which 
can  be  manifested  as  variations  in  the  local  coercive  held  and  local  remanent  polarization.  Other 
variations  can  be  produced  by  stress  nonhomogeneities  and  variable  effective  fields. 

To  incorporate  these  effects  en  route  to  constructing  a  macroscopic  model,  we  consider  the  coercive 
held  Ec  to  be  a  manifestation  of  an  underlying  density  u\  (Ec)  rather  than  fixed  values  which  is  typically 
assumed  for  single  crystals  having  a  uniform  lattice  structure.  To  create  a  macroscopic  model  for 
the  polarization,  we  also  consider  the  variation  of  effective  fields  in  the  material.  As  detailed  in 
[6,  8],  an  applied  held  £  in  a  ferroelectric  material  is  augmented  by  an  interaction  held  Ej  generated 
by  neighboring  dipoles  which  produce  nonhomogeneous  effective  fields  in  the  material.  This,  along 
with  various  other  processes,  produces  variations  in  the  applied  held  that  can  signihcantly  alter  the 
resulting  polarization.  To  incorporate  these  variations,  we  consider  the  effective  held  Ee  =  E  +  Ej  to 
be  distributed  about  the  applied  held  E  with  an  underlying  density  for  the  interaction  held  Ej  which 
we  denote  by  h/2{ Ej ).  The  introduction  of  variations  in  the  effective  held  produces  domain  switching 
in  advance  of  the  remanence  point  in  accordance  with  observations  from  experimental  data. 

The  complete  macroscopic  polarization  model  for  nonhomogeneous,  polycrystalline  materials  with 
distributed  coercive  and  effective  helds  is 

/*oo  /*oo 

P(E)=  /  /  {P){E-EI,Ec)ul{Ec)v2{EI)dEIdEc  (12) 

J  0  J — 00 

where  (P)  is  the  local  polarization  kernel  given  by  (5)  when  thermal  effects  are  negligible  or  by  (11) 
when  incorporating  thermal  relaxation. 
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5  Model  Validation  and  Properties 


In  this  section,  the  ferroelastic  model  presented  in  Section  3  is  compared  with  experimental  PLZT 
data  reported  in  [2] .  The  ferroelastic  switching  model  employs  the  thermal  activation  model  presented 
in  Section  3.3  to  evolve  the  polarization  in  time.  However,  since  the  measured  values  exhibit  minimal 
relaxation,  the  switching  model  was  computed  with  a  large  7  and  a  small  r  which  resulted  in  negligable 
thermal  activation.  As  reported  in  [2],  the  PLZT  composition  is  also  very  near  the  rhombohedral- 
tetragonal  morphotropic  boundary.  To  characterize  this  nonhomogenous  composition,  we  utilize  the 
macroscopic  polarization  model  of  Section  4. 

For  the  macroscopic  model,  we  make  the  a  priori  assumption  that  the  density  for  the  coercieve 
field  is  given  by 


vi(Ec(a)) 


Cie-E2c(<r)/bl  e-El(a)/bl  >  q 
0  g-BeW/H  <  0 


(13) 


and  the  density  for  the  interaction  field  is  given  by 


v-2 (£/)  =  c2e  Ei/b2. 


(14) 


The  relation  Ec(a)  is  specified  by  the  single  crystal  value  (6)  and  the  variance  of  the  coercive  field 
was  estimated  from  the  variance  in  the  experimentally  measured  values  in  Table  2.  The  variance  of 
the  effective  field  was  estimated  by  the  degree  to  which  switching  occurs  before  remanence  in  the 
PLZT  data.  The  parameters  of  the  densities  17 (Ec)  and  17 (Ej)  used  in  the  model  were  c\  =  27.0, 
b\  =  6.5  x  10“4  (MV/rn)2,  C2  =  6.0  and  62  =  5.0  x  10~3  (MV/rn)2.  To  implement  the  macroscopic 
polarization  given  by  (12),  a  composite  Gaussian  quadrature  was  employed. 

The  slope  of  the  hysteron  in  the  180°  regime  is  given  by  ip1 ,  the  slope  of  the  hysteron  in  the  90° 
regime  is  given  by  r p1  and  the  remanent  polarization  is  defined  by  the  parameter  Pr.  The  remaining 
coefficients  a,  q  and  Y  affect  the  slope  and  intercept  of  the  linear  Ec(a )  relation  (6)  and  these  values 
can  be  ascertained  by  fitting  the  Ec(a)  relation  to  experimental  data  given  in  Table  2.  The  values  of 
Pm ,  7 1 ,  /3  and  A  are  defined  in  Equation  (1).  The  value  of  the  parameters  used  to  model  the  PLZT 
data  are  compiled  in  Table  3. 

The  model  is  employed  with  varying  applied  stresses  as  well  as  an  oscillating  electric  field  to  match 
experimental  conditions;  the  reader  is  referred  to  [2]  for  details  regarding  the  experimental  procedures. 
The  ferroelastic  model  is  compared  to  the  PLZT  data  in  Figures  3-5.  The  model  parameters  were 
chosen  to  optimize  the  fit  of  the  E  —  P  and  E  —  e  data  shown  in  Figures  3  and  4.  The  ferroelastic 
model  characterizes  the  90°  switching  that  occurs  in  the  E  —  P  data  as  a  compressive  stress  is  applied 
to  the  PLZT  sample.  The  model  also  predicts  a  negative  strain  due  to  an  applied  compressive  stress 
as  well  as  the  butterfly  nature  of  the  E  —  e  data. 

It  is  noted  that  a  possible  source  of  error  in  the  model  fit  may  arise  from  the  rhombohedral  nature 
of  the  PLZT  data  since  the  ferroelastic  model  is  derived  in  the  tetragonal  phase. 

The  datafit  of  the  <7  —  e  relation  shown  in  Figure  5  can  be  optimized  by  setting  the  Young’s  modulus 
parameters  to  match  the  slope  of  the  appropriate  part  of  the  a  —  s  curve.  However,  this  results  in  an 


a(MPa) 

PR(C/m2) 

Ec(MV/m) 

Ef  (MV/m) 

0 

0.247 

0.35  ±  0.01 

0.35  ±  0.01 

-6 

0.247 

0.26  ±  0.01 

0.35  ±  0.01 

-10 

0.235 

0.19  ±  0.02 

0.35  ±  0.02 

-15 

0.215 

0.15  ±  0.02 

0.39  ±  0.02 

Table  2:  Experimental  values  for  <7,  Pr,  E®°,  and  Ec  for  rhombohedral  PLZT  from  data  reported  in  [2]. 
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parameter 

value 

parameter 

value 

Pr 

0.24  (C/m2) 

a 

0.01  (MV/rn) 

Pi 

0.239  (C/m2) 

q 

500.0  (MV/rn)2 

P901 

0.015  (C/m2) 

Y 

20000  (MPa) 

V 

20.0 

V  2 

20.0 

Table  3:  Parameters  for  the  1-D  piecewise  ferroelastic  switching  model. 


underprediction  of  the  strains  in  the  E  —  e  curve.  Simultaneously  optimizing  both  the  a  —  e  and  E  —  e 
relations  may  be  possible  by  using  higher-order  terms  in  the  electromechanical  coupling  energy  given 
by  (2). 

The  PLZT  data  in  [2]  illustrates  the  stress-dependence  of  the  remanence  and  saturation  polariza¬ 
tion.  For  applied  compressive  stresses  greater  than  15  MPa,  the  remanence  and  saturation  polarization 
significantly  decrease.  Since  the  remanence  polarization  is  a  fixed  model  parameter,  the  present  for¬ 
mulation  of  the  ferroelastic  model  should  be  limited  to  moderate  stress  regimes. 


Figure  3:  Macroscopic  E  —  P  relation  for  varying  applied  stresses.  Homogenized  1-D  ferroelastic  model 

( - )  and  experimental  PLZT  data  (•  •  • ):  (a)  <73  =  0  MPa,  (b)  <73  =  —6  MPa,  (c)  <73  =  —10  MPa, 

and  (d)  03  =  —15  MPa. 


o3  =  0  MPa  o3  =  -6  MPa 


-0.6  -0.4  -0.2  0  0.2  0.4  0.6  -0.6  -0.4  -0.2  0  0.2  0.4  0.6 


Eg  (MV/m)  Eg  (MV/m) 

(c)  (d) 

Figure  4:  Macroscopic  E  —  e  relation  for  varying  applied  stresses.  Homogenized  1-D  ferroelastic  model 

( - )  and  experimental  PLZT  data  (•  •  • ):  (a)  <73  =  0  MPa,  (b)  <73  =  —6  MPa,  (c)  (73  =  —10  MPa, 

and  (d)  <73  =  —15  MPa. 

6  Concluding  Remarks 

This  paper  summaries  a  model  characterizing  stress-induced  90°  and  180°  switching  inherent  to  lead 
zirconate-titanate  in  a  manner  suitable  for  subsequent  transducer  and  control  design.  Helmholtz  and 
Gibbs  energy  functionals  were  used  to  characterize  the  electromechanical  behavior  of  homogeneous, 
single  crystal  compounds.  The  construction  is  similar  to  that  employed  for  SMA  [3,  5,  6,  7]  in  the 
sense  that  it  is  a  1-D  potential  with  three  wells  corresponding  to  the  ±180°  and  90°  equilibria.  The 
construction  of  this  functional  is  phenomenological  but  the  resulting  decrease  in  dimension  signifi¬ 
cantly  diminishes  implementation  time.  The  functionals  are  directly  minimized  to  provide  kernels  for 
characterization  in  the  absence  of  thermal  relaxation  or  balanced  with  the  relative  thermal  energy 
through  Boltzmann  principles  to  incorporate  relaxation  phenomena.  The  effects  of  material  nonho- 
mogeneities  and  poly  crystallinity  are  incorporated  by  assuming  that  properties  such  as  local  coercive 
and  interaction  fields  are  manifestations  of  underlying  distributions  rather  than  constants.  Stochastic 
homogenization  in  this  manner  provides  a  low  order  macroscopic  model  suitable  for  transducer  and 
control  design. 
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Figure  5:  (a)  Macroscopic  a  —  P  relation  and  (b)  macroscopic  a  —  e  relation.  Homogenized  1-D 
ferroelastic  model  ( - )  and  experimental  PLZT  data  (•••). 
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